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ABSTRACT 

Recent observations show a large concentration of galaxies at high redshift. At first 
sight strong clustering of galaxies at high redshifts seems to be in contradiction with 
the models of structure formation. In this paper we show that such structures are a 
manifestation of the strong clustering of rare peaks in the density field. We compute 
the frequency of occurrence of such large concentrations of galaxies in some models of 
structure formation. 
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1 INTRODUCTION 

Observations of galaxies at hig h redshifts, around 2 = 3, 
in a contiguous region of sky (steidel et al 1997) show a 
marked concentration in a region of width Az = 0.04, or 
Av = 3000km s _1 . The angular size of the region probed is 
9' x 18' and within this region the galaxies contributing to 
the peak are distributed randomly in about half the area. A 
quasar is also present in this wall like structure. The clus- 
tering in redshift space is estimated to have a confidence 
level of 99.8%. Steidel at al. (1997) point out that cluster- 
ing at such large scales, in most models of galaxy formation, 
requires high bias to explain the observed structure. 

First galaxies to form in the universe correspond to the 
deepest potential wells or the highest peaks in the initial 
density distribution. It is well known that these ra re peaks 
cluster mor e strongly than, say, the typical peaks (Bardeen 



et. al. 1986). Therefore, at early times, when only rare peaks 



have collapsed into structures like galaxies, we expect these 
objects to show significant clustering. The epoch when typi- 
cal halos of a given mass scale M collaps e is characterised by 
o~(M, z) ~ 5 C (Press and Schechter 1975). It is customary to 
use S c = 1.69 as a spherical perturbation virialise s when the 
linearly ex trapolated density contrast equals S c (Gunn and 
Gott 1973). If we define a quantity v(M,z) = S c /a(M,z), 



then we can say that typical halos of mass M collapse when 
v(M) ~ 1. For v 3> 1, only the rare peaks of the given mass 
scale have collapsed and these tend to cluster very strongly. 
In this paper, we study clustering of halos at z ~ 3 and com- 
pare it with the constraint set by the observation of Steidel 
et al. (1997) for some models of structure formation. 

We assume, that the galaxy distribution and the halo 
distribution are identical. This may not be true, especially at 
small scales where gastrophysical processes play an impor- 
tant role. However, including these processes requires more 



detailed modelling which, in turn, requires more assump- 
tions. Therefore, we will restrict the present study to the 
study of halo distribution at high redshifts. 
We address the following questions here: 

• Is the clustering of dark matters halos in models of 
structure formation sufficiently strong to explain the ob- 
served concentration of galaxies? 

• What is the frequency with which we may expect to 
see a 3<r excess for the number of halos if we simulate the 
observations of galaxies using numerical simulations? 



2 CLUSTERING OF HALOS 

We can study the clustering of halos using numerical simula- 
tions. For this particular study we ch oose the SCDM model 
with the sha pe parameter T = 0.5 (Efstathiou, Bond and 
White 1992). We normalise the power spectrum so that 
as — o"(8h _1 Mpc, z = 0) = 0.6. Most models have more 



power at the relevant scales as compared to this, and if this 
model provide enough clustering at z = 3, most other mod- 
els can do that too. 

The observed number density of Lyman break galax- 
ies is 2.9 x 10 -3 h 3 Mpc~ 3 in the Einstein-de Sitter Universe. 
This im plies a halo mass of about 5 x 10 12 M Q for the SCDM 
model ( |Press and Schechter 1975 ) . Using this mass for halos 
is same as assuming that all halos host Lyman break galax- 
ies. In order to relax this assumption, we choose a lower 
halo mass, M^aio > 7x 10 11 Mq. Thus less than one in ten of 
these halos host Lyman break galaxies. In absence of any al- 
ternative, we assume that galaxies are distributed randomly 
amongst these halos and hence have the same clustering 
properties. 

We quantify the clustering of halos with the averaged 
two point correlation function. This is defined as 
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Figure 1. This figure shows the averaged two point correlation 
function for halos of mass M^alo > 7 x 10 Mq at z = 3. The x- 
axis shows the comoving scale in units of h — iJVlpc. The amplitude 
of the halo correlation function is comparable to the galaxy corre- 
lation function at z = 0, i.e. it is much higher than the amplitude 
of mass correlation function at z = 3. 



x ( > {x)dx - 



3J 3 (r) 



(1) 



where £(a;) is the two point correlation function. 

Halos are identified using the Friends-of Friends (FOF) 
algorithm with a linking length of 0.2 (over-density of 60). 
Fig.l shows the averaged correlation function as a function 
of scale for these halos at redshift z = 3. If the number 
density of these halos is n, then the number of halos within 
distance r from a typical halos is 47rr 3 n(l + ^(r))/3. Thus we 
expect to see twice the average number in a sphere of radius 
8h~ Mpc around a typical halo. Some halos will have more 
neighbours than a typical halo. Therefore, an excess of a 
factor four in one bin is not very surprising. 

The large amplitude of correlation function is not com- 
pletely surprising. Analytical models for evolution of bias for 
halos also predict that bias increases with increasing redshift 
(For example, see Matarrese et al. 1997). These predictions 
suggest that the amplitude of halo correlation function de- 
creases very slowly with increasing redshift. However, most 
of these models do not predict the halo correlation function 
at z — 3 t o be as stron g as the present day galaxy correlation 
function flBagla 1997|). 



3 SYNTHETIC OBSERVATIONS 

In this section, we will simulate observations using N-Body 
simulations. We will restrict our study to the CDM class 
of models. In particular, we will study clustering in the 
standard CDM model (T = 0.5) and a variant of SCDM 
that reproduces the observed correlation function out to 



large scales (F = 0.3). These models are normalised so that 
as = a(8h~ 1 Mpc, z = 0) = 0.6. To describe differences in- 
duced by a higher normalisation we will use the same models 
with as — 1. The background cosmology is assumed to be 
S7 = 1, Ho = 50km Mpc - s _1 . We will specify all scales in 
the comoving coordinates. 

The observations indicate a very large scale for the 
concentration of galaxies. In the Einstein deSitter model 
Az = 0.04 translates to 15ft _1 Mpc. The observed angular 
extent of the structure in question is greater than or equal 
to 8' x 11'. (One arc minute corresponds to 0.87/i _1 Mpc at 
z — 3.) To study clustering on such large scales we need 
to simulate a very large volume and we chose to work with 
a simulation box of size 166/i _1 Mpc. All simulations were 
done using 128 3 particles so that mass of each N-Body par- 
ticle was approximately 10 12 Mq. All N-Body simulations 
used here were done using a Particle-Mesh code. 

To simulate "observations" of high redshift galaxies we 
use the following method. 

• As galaxies form in regions with high density, we begin 
by isolating such regions in the simulation volume. 

• We project the particles in the selected regions in red- 
shift space by assuming that one of the axes is aligned along 
the line of sight. Velocity along the line of sight and the 
Hubble redshift combine to give the total effective redshift 
of each particle. 



1 + Ztot = (1 + Z hub ) (1 + Z pe 



(2) 



Here z pec = v p /c with c the speed of light and v p the com- 
ponent of the peculiar velocity of the particle along the line 
of sight. 

• We then view the selected particles around a large num- 
ber of lines of sight. For each line of sight all particles within 
a square region of 10' x 10' are included in the field of view. 
The redshift distribution of these particles is then analysed 
for a large number of fields of view. 

• We smooth the distribution in redshift by using a top 
hat window of width Az — 0.04. We search for peaks in the 
smoothed distribution. As the correspondence between par- 
ticles and galaxies is not very clear, we will measure peaks 
from the average amplitude in units of the standard devia- 
tion a. 

We are interested in clustering at large scales and there- 
fore we have to use a large simulation volume. This implies 
that the mass of individual particles in the simulation is 
considerably larger than the mass associated with the galax- 
ies we are studying. This clearly makes it difficult to iden- 
tify galaxies directly and we can only try to isolate regions 
where such galaxies could have formed. Therefore, we work 
with three density thresholds: (1) g c = 0, i.e. no cutoff; (2) 
q c — 2gt, i.e. we will select regions with S > 1. Here Qt is the 
background density and S is the density contrast. And, (3) 
q c = 4gt, or S > 3. As we are using very low thresholds, the 
frequency of occurrence of large concentrations of galaxies 
computed here is essentially a lower limit. 

Our aim here is to search for clustering at large scales in 
redshift and the width of the structure we hope to reproduce 
is Az = 0.04. However, it is useful to start with smaller bins 
and average over the neighbouring bins with a top hat filter. 
This ensures that all high peaks are picked out. 

We group the particles in each field of view in bins of 
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Figure 2. This figure shows the dense structures selected using 
a threshold g c = ^Qb tn real as well as redshift space. The model 
used in these simulations is CDM with T = 0.5, normalised to 
erg = 0.6. The figures show projected distribution of a slice of 
thickness 10'. The top panel shows the particles in real space 
and for comparison with the second panel and observations, we 
have used the cosmological redshift and angle cj> as coordinates. 
The lower panel shows the same set of particles in redshift space. 
Here 1 + z to t = (1 + z hub )(l + z pec ), where z pec = v p /c with c as 
the speed of light and v p the component of peculiar velocity along 
the line of sight. It is obvious from this figure that redshift space 
distortions operate on much smaller scales than Az = 0.04. It is 
also clear that there are clearly some structures with an angular 
extent greater than 10', size of the observed concentration. 



500km s _1 . We then combined six bins around each of the 
smaller bins to get a distribution with Az — 0.04 sampled 
at every 500km s _1 . We scan this distribution for locating 
maxima and in order to avoid over-counting high peaks, we 
ensure that a given bin has the largest number of particles 
as compared to bins within Az — 0.02 on either side. The 
amplitude of peaks is measured from the average in units 
of standard deviation - the average and standard deviation 
are obtained from a large number of "fields of view" through 
the given simulation box. 



4 RESULTS 

We have plotted a slice from an N-Body simulation in fig. 2. 
This figure shows the projection of a slice from a CDM sim- 
ulation with r = 0.5 and as = 0.6. The particles shown 
here are located in regions with 8 > 3. The projection is 
shown in both real and redshift space. For easy comparison 
with the relevant observations, we have used angular size 
and redshift as the coordinates. The thickness of this slice is 
10'. The top panel shows the distribution of particles in real 
space and the lower panel shows the same set of particles in 
redshift space. It is clear that redshift space distortion tends 
to squeeze many structures into a thin sheet perpendicular 
to the line of sight. However, the scale where this effect is 
important is much smaller than Az = 0.04. The sizes of 
largest structures - angular sizes as well as the extent in 
redshift - and the effects of redshift space distortion can be 
seen clearly. Redshift space distortions play a very impor- 
tant role at small scales (McGill 1990a, 1990b) but for the 
scales that we are interested in, i.e. Az — 0.04, these effects 
are smoothed out. At these scales, the frequency with which 
such structures appear is identical in the real and redshift 
space. 

Fig. 3 shows the redshift distribution in two fields of view 
obtained using the method outlined in §3. These particular 
examples were chosen to exhibit the variety of distributions 
and peaks seen in such systems. In one case we see a semi- 
periodic set of peaks whereas in the other case we see only 
one isolated peak. 

Table 1 lists the frequency of occurrence of Na peaks 
above average in the redshift space distribution. Here we 
have extrapolated from the range of redshift covered by the 
simulation to that used in Steidel et al. (1997) to facilitate 
comparison. We can summarise the conclusions as follows: 

• Models with same normalisation predict similar fre- 
quency of occurrence for peaks with Na more particles than 
average. This is to be expected as both models have similar 
power at the relevant scale. 

• Models with higher normalisation predict a higher fre- 
quency of occurrence for the unbiased case. 

• Models with higher normalisation predict a lower fre- 
quency of occurrence for 5a peaks above average for the high 
density threshold. This can be understood if the statistical 
bias depends only on v defined in the introduction. This is 
small for models with the higher normalisation, and hence 
we do not see the effect of strong clustering of rare peaks. 
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Figure 3. Distribution of objects in redshift space. These pan- 
els show the redshift distribution in two 10' X 10' fields of view 
through a simulation volume. These are from the simulation used 
in fig. 2. These frames show a few high peaks in the distribution. 
The y axis is the height of a peak measured from the average in 
units of standard deviation a. The curve shows an average over 
Az = 0.04 sampled at intervals of Az = 0.0067. 



5 DISCUSSION 

We have shown that the observations of strong clustering of 
galaxies at high redshifts are not in conflict with the pop- 
ular models of structure formation. We have demonstrated 
this in two different ways: first by showing that the bias 
can indeed be very large at early times. We have also gener- 
ated synthetic observations in a few models and have shown 
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Table 1 . This table lists the frequency of Na peaks above average 
in the models of structure formation studied here. The first two 
columns refer to the parameters of the model, the shape parame- 
ter r and normalisation erg. The third column lists the threshold 
in density contrast used to select particles. The last three columns 
list the frequency for finding 3, 4 and 5 a peaks above average in 
fields of view of 10' X 10' and redshift range 2 < z < 3.5. 



that the frequency of occurrence of large concentrations of 
galaxies is compatible with observations. 

Fig. 2 suggests that still larger structures maybe found 
in future searches. The distribution of sizes of largest struc- 
tures in surveys at high redshift - both angular size and the 
extent in redshift - may be used to discriminate between 
different models. 

Strong clustering of galaxies at high redshifts implies 
that the galaxy clustering evolves in a very different man- 
ner as compared to clustering in the underlying mass distri- 
bution. A detailed study of the evolution of halo clustering 
and its imp lications for galaxy clustering has been presented 
elsewhere (Bagla 1997). 
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